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Abstract 

We use CALIOP nighttime measurements of lidar backscatter, color and depolarization ratios 
during the summer of 2007 to study transatlantic dust properties downwind of Saharan 
sources, and to examine the interaction of clouds and dust. Our analysis suggests that (1) 
while lidar backscatter doesn’t change much with altitude in the Saharan Air Layer (SAL), 
depolarization and color ratios both increase with altitude in the SAL; (2) lidar backscatter 
and color ratio increase as dust is transported westward in the SAL; (3) the vertical lapse rate 
of dust depolarization ratio, introduced here, increases within SAL as plumes move westward; 
(4) nearby clouds barely affect the backscatter and color ratio of dust volumes within SAL but 
not so below SAL. Moreover, the presence of nearby clouds tends to decrease the 
depolarization of dust volumes within SAL. Finally, (5) the odds of CALIOP finding dust 
below SAL next to clouds are about 2/3 of those far away from clouds. This feature, together 
with an apparent increase in depolarization ratio near clouds, indicates that particles in some 
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dusty volumes lose asphericity in the humid air near clouds, and cannot be identified by 
CALIPSO as dust. 

1 Introduction 

Atmospheric mineral dust particles have significant effects on the climate and the 
environment. Despite notable advances in modeling and satellite and ground-based 
measurements, dust remains the dominant factor in the uncertainty of aerosol radiative forcing 
(IPCC, 2001, 2007). Dust emitted from dry areas of Africa is transported over the North 
Atlantic Ocean to coastal areas of America with a peak of deposition during summer months 
(e.g., Prospero and Carlson, 1972; Mattsson and Nihlen, 1996; Prospero and Lamb, 2003; 
Torres et al., 2002; Kaufman et al., 2005). 

The influence of dust on a radiative budget depends on its ability for absorbing and scattering 
solar and IR radiation. Dust optical properties are determined by the refractive index (i.e., 
chemical composition) (e.g, Sokolik and Toon, 1999; Wang et al., 2002; Lafon et al., 2006; 
Kahnert et al., 2007; Kandler et al., 2007, 2009; Osborne et al., 2008; Petzold et al., 2009), 
size and shape of dust particles (Kalashnikova and Sokolik, 2002; Dubovik et al., 2006; 
Nousiainen, 2009). Dust hygroscopicity describes the particles’ ability for taking up water 
from humid air, and is rooted in the physical-chemical properties of dust components. Dust is 
mostly composed of water-insoluble minerals and shows nearly complete hydrophobicity or 
at least poor hygroscopicity (e.g. Twomey, 1977; Li-Jones et al., 1998). However, after being 
lifted in the air and mixed/coated with water-soluble materials such as sea-salt, sulfate, or 
nitrate by atmospheric processing (Levin et al., 1996; Yin et al., 2002), dust hygroscopicity 
can increase, which can cause changes in optical properties. The increase in hygroscopicity 
may also transforms dust into effective cloud condensation nuclei (CCN) or ice nuclei (IN) 
(Johnson, 1982; Wurzler et al., 2000; DeMott et al., 2003; Sassen et al., 2003; Lohmann and 
Feichter, 2005; Twohy et al., 2009a) and thereby affect the formation and distribution of 
clouds and precipitation (e. g. Kelly et al., 2007), thus also altering the radiative impact of 
clouds. 

The evolving shape-dependent optical and hygroscopic properties of dust pose a question: 
How does the hygroscopicity of dust affect its optical properties near clouds? This question is 
pertinent as recent studies have shown that optical properties of clear sky aerosols different in 
the vicinity of clouds from those far away from clouds (e.g., Clarke et al. 2002; Twohy, et al, 
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2002, 2009b; Koren et al., 2007; Su et al., 2008; Redemann et al., 2009; Tackett and 
Girolamo, 2009; Varnai and Marshak, 2011). The answer to this question is likely not only to 
improve our understanding of dust-cloud interactions but also yield better estimates of direct 
radiative forcing. To that end, we present analysis of dust properties over the North Atlantic 
Ocean — including near-cloud behavior — based on Cloud-Aerosol Lidar with Orthogonal 
Polarization (CALIOP) data (Winker, et al. 2003). 

CALIOP is a space based lidar system onboard the Cloud Aerosol Lidar and Infrared 
Pathfinder Satellite Observations (CALIPSO) satellite launched in 2006. CALIOP data offers 
many advantages for this study. First, since CALIOP uses a laser with a small footprint (~90 
m in diameter on the ground), its aerosol data is not affected by the 3D radiative 
enhancements of nearby clouds, which cause complications for instruments observing 
reflected sunlight (Wen et al., 2007; Marshak et al., 2008; Varnai and Marshak, 2011). 
Second, CALIOP provides backscatter depolarization information at 532 nm, which allows 
one to distinguish (typically) non-spherical dust from (typically) spherical droplets (Sassen, 
2000; Murayama et al., 2001 Vaughan et al., 2004). Third, CALIOP’s high spatial resolution 
(30 m vertically and 333 m horizontally) is well-suited for studying cloud-dust interactions 
that have typical scales of several kilometers (e. g., Koren et al., 2007; Varnai and Marshak, 
2011). 

2 Data and methodology 

This study uses a month-long (06/07/2007-07/07/2007) dataset of nighttime CALIOP Version 

3 data over the North Atlantic Ocean (0°-45 North, 0°-90° West). We note that since 
CALIPSO orbits are repeated in a 16-day cycle, our month-long dataset covers almost two 
orbital cycles, with the longitudes of closest orbits being 1.55° apart. Since in summer there 
are usually multiple outbreaks per month (Huang et al., 2010), this month-long dataset is 
sufficient to observe the basic features of dust outbreaks. 

CALIOP measures the total backscatter of its laser pulses at 532 nm and 1064 nm 
wavelengths, and the perpendicularly polarized backscatter at 532nm. 

This paper characterizes dust particles using their attenuated backscatter coefficient (|3'), 
attenuated color ratio (%'), and depolarization ratio (5') values. Unless specified otherwise, (’>' 
is the median of the vertically averaged 532 nm attenuated backscatter coefficients within 
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dust layers identified by the 5 km-resolution aerosol product. %' is the ratio between values 
at 1064 nm and 532 nm wavelengths, which tends to increase with the size of particles; this is 
especially true for optically small spherical particles. 5' is the ratio between the perpendicular 
and parallel components of P'532 , and is greater for non-spherical particles. We note that 
analyzing measured optical quantities helps avoid potential artifacts that could be introduced 
by uncertainties in retrieving physical particle properties. However, we still use the results of 
operational CALIOP data processing for identifying cloud and dust layers. The operational 
algorithm identifies cloud and dust layers in three steps. 

First, it identifies particle layers based on the observed 532 nm backscatter values (Winker et 
al., 2009; Vaughan et al., 2009). 

Second, it determines whether a detected layer is a cloud or aerosol layer based on its latitude, 
altitude, 532 nm backscatter, color ratio and depolarization ratio (Liu et al., 2004, 2009). The 
most obvious clouds are identified at 333 m resolution, while the more ambiguous cases are 
decided at a coarser (1 km or 5 km) resolution. Our study considers a location cloudy if the 1 
km resolution Level 2 CALIOP cloud product indicates the presence of clouds. In order to 
reduce the impact of misclassifications between clouds and aerosols, this paper examines 
aerosol layers only if the Cloud-Aerosol Discrimination (CAD) product (Liu et al., 2004, 
2009) — based on probability distribution functions obtained from expert classifications for 
sample orbits — indicates that the likelihood of aerosol exceeds 70%. As additional precaution, 
the statistical analysis in this paper examines median values instead of mean values, because 
medians are less sensitive to any outlying data points influenced by undetected cloud 
particles. The uncertainty of median values is estimated using the bootstrapping algorithm 
(Efron and Gong, 1983). 

Third, the operational CALIOP data processing identifies dust layers as aerosol layers with 
high depolarization ratio values (Omar et al, 2009). Because depolarization depends on 
particle shape, it is well suited for separating typically non-spherical dust particles (5' > 0.2) 
from usually spherical non-dust aerosols (5' < 0.075) over ocean. We note that this paper 
considers dust-containing layers identified as either “dust’ or “polluted dust” in the 5 km- 
resolution Level 2 aerosol product. (Polluted dust is dust mixed with biomass burning 
aerosols or polluted marine aerosols, with a depolarization ratio between those of dust and 
non-dust aerosols, 0.075 and 0.2 (Omar et al., 2009)). 
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To discern changes in dust properties during transatlantic transport, we examine dust behavior 
in the three regions shown in Figure 1: East (E) (0°-30° W), Middle (M) (30°-60° W) and 
West (W) (60°-90° W). These three regions lie at different distances from the African dust 
sources, and cover most of the dust paths from Africa to America during the summer of 2007. 

3 Results and discussion 

The spatial and optical characteristics of African dust vary during the transatlantic journey (e. 
g., Liu et al., 2008; Huang et al., 2010). In this section we examine the variations in three 
steps. First, the overall statistics of dust properties in the three geographic regions are 
compared. This part focuses on the vertical distribution of dust samples in the CALIOP 5 km 
resolution aerosol product, and on the vertical distributions of attenuated backscatter 
coefficient, color ratio and depolarization ratio. We then analyze the relationships between 
dust properties and cloud coverage in the three regions. Finally, we discuss the systematic 
changes in dust properties that occur near clouds. 

3.1 Dust properties in the three regions 

3.1.1 Vertical distribution of dust 

Figure 2a shows the vertical distribution of dust samples in our three regions, normalized by 
the total number of dust samples within each region. The number of dust samples is defined 
as the number of 5 km long 270 m high volumes that, according to the CALIOP aerosol 
product, contain dust and have a CAD value between -70 and -100. In the Eastern (E) region, 
more than 80% of dust is between 1.5 km and 5.5 km altitude, with the peak probability 
around 3.5 km. This elevated dust distribution is a typical result of two confining inversions 
below and above the SAL (Carlson and Prospero, 1972). The dust remains elevated in the 
middle (M) region as well, although the mean elevation descends about 0.5 km. The similar 
patterns in the E and M regions indicate that the driving forces keeping the dust at high 
altitudes in region E persist in region M as well. If a 3-day average transport time from 
Region E to M is assumed, the descending velocity from center of E to center of M is 
estimated around 1.7 mm/s, which is consistent with the typical SAL average descending 
velocity of 1-2 mm/s (Carlson and Prospero, 1972). Finally, in the West (W) region the 
chances of finding dust decrease steadily with altitude, and dust is rarely found above 5 km. 
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This dramatic change in the vertical distribution implies that the meteorological conditions in 
region W are different from those sustaining the elevated profdes in regions E and M. 
Moreover, the vertical distribution of dust in region W indicates that (dry and wet) dust 
sedimentation has the strongest impact over region W. 

3.1.2 Attenuated backscatter coefficient, color ratio, and depolarization 
ratio of dust 

This section examines dust optical properties in the three regions. Since dust above 5 km is 
rare, the analysis of dust optical properties will be limited to altitudes below 5 km. The 
analysis uses the Level 2 CALIOP aerosol layer product, which provides averaged (3', 7 ', and 
5' values for each dust layer. These layer-average values are assigned to all altitude bins 
within a dust layer when creating Figure 2 (and 4). 

The results in Figure 2 show that dust properties vary with altitude differently within SAL 
(from 1.5 km to 5 km in altitude of regions E and M), below the SAL (below 1.5 km in 
regions E and M), and in region W. 

In the SAL, median P '532 values are nearly constant with altitude, but the medians of 7 ' and 5' 
increase with altitude. Since the typically large and non-spherical dust particles imply large 7 ' 
and 5' values, the observed dusty volume behavior suggests that the concentration of dust 
increases with altitude and/or the concentration of non-dust marine aerosols in dusty volumes 
decreases with altitude. 

Below the SAL in the E and M regions, the median 5' increases with altitude. This is the 
result of dust mixing with non-dust marine aerosols in the moist air confined between the 
marine surface and the inversion created by the dry and warm SAL aloft. Since the 
concentration of wet marine aerosols is much higher below than inside the SAL, backscatter 
from these aerosols contributes significantly to the lidar signals and reduce the depolarization 
ratios of dusty volumes below the SAL. As discussed in the following paragraph, the mixing 
of dust and marine aerosols in region W could be the major reason behind higher backscatter 
and lower depolarization in region W than in E or M. 

In region W, the median of P '532 decreases with altitude and the medians of 7 ' and 5' increase 
with altitude up to at least 4-5 km. This behavior arises from the dissipation of the SAL in the 
region. As the SAL extends westward, its temperature drops and convection can bring in 
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moist air from below. Eventually the inversion confining the moist marine air at low altitudes 
breaks up in region W, which allows the moist marine air to reach much higher altitudes. 
Even so, the concentration of wetter and larger marine aerosols tends to decrease with 
altitude, and so their contribution to the lidar backscatter of dusty volumes tends to be smaller 
at higher altitudes. This results in higher median x' and 5' values at higher altitudes even in 
region W. 

A comparison of dust properties in regions E, M, and W also reveals several features related 
to SAL influences at different transport stages. For example, at most altitudes the medians of 
P '532 and x' increase westward, whereas the median of 5' tends to decrease westward. These 
trends arise from the fading of the SAL during the westward transport: As the air moves from 
west Africa to the Caribbean, the SAL temperature keeps decreasing, which allows the 
concentration of moist marine air mixed in from below to keep increasing. The presence of 
more moist marine aerosols at west increases the median P '532 and 7 / of dusty volumes and 
reduces the median 5' value. 

Flowever, some features in Figure 2 cannot be explained by contributions from non-dust 
marine aerosols in dusty sample volumes, and are likely caused by changes in the properties 
of dust particles instead. For example, Figure 2d shows that above 3.5 km, the median 5' is 
larger in region M than in region E. This cannot be explained by mixing from below, because 
the larger contribution of mixed-in marine aerosols in region M would imply smaller (rather 
than larger) 6 ' values. Instead, the observed tendency is likely related to lower fall speed for 
aspherical dust particles: as the more spherical particles fall faster, this leaves an increasingly 
non-spherical dust population at high altitudes as the air moves to region M. The plausibility 
of this scenario is also supported by simulations for highly irregular particles falling slower 
than more spherical ones, because of greater air resistance (Ginoux, 2003). This issue will be 
further discussed in Section 3.3.2. 

We note that the slight increase of x' of dusty volumes from E to M (in Figure 2c) is different 
from the behavior of dust Angstrom exponents retrieved by MISR and MODIS, the latter 
displaying no significant changes during transatlantic transport (Kalashnikova and Kahn, 
2008). This contrast is about vertical resolution: the CALIOP trend of westward increase in 
Figure 2c occurs within SAL but MISR and MODIS integrate the entire atmospheric column 
and hence are affected by particles above and below SAL. 
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3.2 Correlation of dust properties with cloud fraction 

This section examines the relationships between dust properties and cloudiness in the three 
study regions. We characterize cloudiness through the cloud fraction (CF), defined for each 
dust-containing 5 km-size column as the ratio of number of cloudy 0.333 km profiles to the 
total number of 0.333 km profiles in the column. Simply put, if the number of cloudy 0.333 
km profiles is m, the cloud fraction is m/\ 5. We note that in addition to the cloud fraction 
varying between 0 and 1 , the relative location of dust and clouds within 5 km wide columns 
can also vary (Figure 3). We also note that unlike the conventional cloud fraction that is based 
on 2-dimensional (2D) images, our definition here is based on 1 -dimensional (ID) 
measurements along the CALIPSO track. Although off-track clouds may influence dust 
properties along the track even for CF\ D =0 (Varnai and Marshak, 2012), CF id is still a 
generally useful indicator of cloud coverage. 

The results in Figure 4 show that dust properties are closely related to CF in all three regions. 
The main features of the relationship are as follows. 

First, the top row of Figure 4 reveals that a smaller fraction of dust samples occurs under clear 
skies in region M than in region E. This is because the SAL is warmer and drier in the East, 
and so the conditions are less favorable for cloud formation in region E than region M. 

Second, rows 2 and 3 in Figure 4 reveal that within each region, the median values of P '532 
and y' are larger for higher CFs. This feature is likely caused by aerosols getting hydrated and 
swelling in humid regions containing clouds, although undetected cloud particles may also 
contribute. The figure also shows that in regions E and M, the increase in backscatter and 
color ratio is more pronounced below the SAL than inside it. The swelling is greater below 
the SAL than inside it both because clouds and high humidity are more common below the 
SAL, and because hygroscopic marine aerosols are fairly abundant at low altitudes even in 
dust layers, whereas the SAL is dominated by less hygroscopic dust particles. 

Third, within each region, two opposite trends of correlations appear between 8 ' and CF: 
Inside the SAL, the median 8 ' of dust is always larger in clear sky than in cloudy skies; 
whereas below the SAL, the median 8 ' of dust is always smaller in clear sky than in cloudy 
skies. The domains of these opposite behaviors can be separated in the fourth row of Figure 4 
roughly at the crossing point of the red curve (CF= 0) and the blue curve (0<CF<0.6). We 
note that these crossing points are approximately at the altitude of the bottom of the SAL. The 
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opposite trends inside and below the SAL clearly indicate a different dust volume 
depolarization response to increased humidity. The possible mechanisms affecting the 
apparent depolarization ratio of dusty volumes below the SAL will be discussed in Section 
3.4. 

3.3 Features of dust volumes in the SAL under clear skies 

Unlike the dusty volumes below the SAL, where the dust is mixed with humidified non-dust 
aerosols, the SAL is dominated by dust particles. This subsection examines several features of 
dust volumes inside the SAL in regions E and M. To reduce the effects of clouds, dust 
volumes are limited to only those under clear skies. In addition, region W is excluded because 
of its low number of dust samples inside the SAL (Figure 2a). 

3.3.1 Relationship of depolarization ratio with color ratio and backscatter 

Figures 2 and 4 show that while backscatter is fairly uniform vertically, both color ratio and 
depolarization ratio increase markedly with altitude inside the SAL of regions E and M. These 
coinciding increases suggest systematic relationship between depolarization ratio on the one 
hand and color ratio and backscatter on the other (Figure 5). Figure 5a reveals a positive 
relationship between the depolarization ratio and color ratio of dust volumes in the SAL of 
regions E and M. In addition, the dynamic ranges of depolarization ratio and color ratio are 
much wider in region M than in region E. Figure 5b shows a similar relationship between 
depolarization ratio and backscatter, though with a much smaller dynamic range for 
backscatter in region E. Figures 5 c and 5d confirm that similar relationships are valid for a 
different dataset (covering 05/25/2008-06/25/2008) as well. 

The relationships in Figure 5 and the similarity of results from the two independent datasets 
can be attributed to the steady altitude-dependence of backscatter, color ratio and 
depolarization ratio inside the SAL. As discussed in Section 3.1.2, these altitude dependences 
are likely caused by two mechanisms: (?) a decrease with altitude in the concentration of non- 
dust particles mixed in from below, and (/'/) different fall speeds vertically separating the 
relatively more spherical dust particles from the least spherical ones. This latter mechanism 
dominates and is further explored in Section 3.3.2. 

In addition, we note that the depolarization ratio is not only a function of the aspect ratio of 
particles but is also related to their size (Mishchenko and Hovenier, 1995). Therefore the 
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relationship between 5' and in Figures 5a and 5c may also indicate that backscatter from 
larger dust particles can be more depolarized. 

3.3.2 Relationship between the vertical increase in depolarization ratio 
and longitude 

As indicated in Figure 2, the depolarization ratio not only increases with altitude in the SAL, 
but also has a larger increase rate in region M than E. As mentioned above, the increase may 
come from more spherical and less spherical dust particles getting vertically separated 
because of their different sedimentation speeds. The upward increase in 5' could then be 
stronger in region M simply because the sedimentation process has more time to work by the 
time the dust reaches region M. 

Numerous studies have demonstrated that the fall speed of atmospheric particles is related to 
their shape (e.g., Ginoux, 2003). Since particles with irregular shapes have greater cross- 
sectional areas and drag-coefficients, they experience stronger drag force in the air — which 
implies that more irregular particles fall slower than more spherical ones. Note that a sphere is 
the most compact object (least surface area for a given volume) and it experiences least drag 
for a given mass. (Here we assume that dust particle shape does not change systematically 
with particle size, which also greatly impacts fall speed.) As a result, shape-induced vertical 
separation will ensue as dust is advected westward, with irregular particles increasingly 
predominant in the upper portions of SAL. At a constant altitude, this stratification is 
expected to widen the dynamic range of depolarization ratios with downstream distance from 
the dust source. This is indeed the case. 

To that end, we divide regions E and M into sub-regions covering 10° wide longitude bands, 
and examine the average difference between the depolarization ratios at 3 and 4 km altitudes 
for each region. As shown in Figure 6 , the average difference 5 ' 4 km - 6 '3 km keeps increasing 
with the distance from the west coast of Africa. This result implies that the observed change 
in volume depolarization ratio within SAL is most likely caused by the greater drag of 
aspherical dust particles. 

3.4 Dust volume properties near clouds 

Relative humidity usually increases as clouds are approached and this causes nearby aerosols 
to swell and get hydrated (acquire thin film of water) or even activated as haze (e.g., Twohy et 
al., 2009b). Observing changes of dust characteristics near clouds can help improve our 
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understanding of the effect of high relative humidity and clouds on dust particles. Figure 4 
has shown that the backscatter and color ratio of dust volumes increase with cloud fraction 
both in and below the SAL, whereas the depolarization ratio changes with cloud fraction 
differently in and below the SAL. This finding indicates that dust properties in and below the 
SAL are different. This section further examines the near cloud behaviors of dust volumes in 
and below the SAL. We note that although the base altitude of the SAL may vary during 
westward transport, as shown in Figures 2 and 4, this analysis uses constant separation 
altitude of 2 km for convenience. The analysis uses CALIOP Level 1 data to examine changes 
in backscatter, color ratio, and depolarization as a function of distance to clouds at a 
resolution of 0.333 km. In this analysis a 0.333 km resolution clear sky profile is considered a 
dusty profile if it is included in one or multiple 5 km resolution dust layer(s). 

Figure 7 illustrates the behavior of dust as a function of distance to clouds. The orange curve 
corresponds to all aerosol samples while the black and green ones to high (in the SAL) and 
low (below the SAL) dust, respectively. Figure 7a shows that the fraction of dust profiles over 
all detected aerosol profiles decreases dramatically near clouds for dust at altitudes below the 
SAL, but remains relatively stable for those in the SAL. The stable behavior in the SAL can 
be explained by the fact that low clouds confined to the boundary layer by the inversion at the 
base of the SAL have little impact on humidity inside the SAL. The near-cloud drop in the 
fraction of dust profiles below the SAL may results from several factors. First, the chances of 
wet removal are higher near clouds, and this can lower the fraction of dusty profiles. Second, 
swelling in the humid air near clouds makes particles more spherical (especially if water- 
soluble particles pollute dust crystals), resulting in the reduction of depolarization ratio thus 
the dust signature (Omar et al., 2009); consequently, some dust-containing profiles are 
(mis)classified as non-dust aerosol. Assuming that below 2 km altitude, the fraction of dust 
profdes is constant beyond 5 km from clouds (as most humidity changes occur within 5 km 
from clouds), the roughly 2/3 drop in the fraction of dust profiles near clouds implies that at 
least 2/3 of dust profiles in the MBL are polluted and hygroscopic. This results in CALIOP 
missing dust in about 1/3 of dust profiles that occur within 5 km from clouds. 

Figures 7b-7d also show that backscatter, color ratio, and depolarization ratio all increase near 
clouds for dust layers below 2 km, but they remain fairly constant for dust layers above 2 km. 
The stable behavior in the SAL occurs because most clouds are below the SAL and have little 
impact on dust in the SAL. In addition, the dust population in the SAL is dominated by 
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hydrophobic particles. For dust volumes below 2 km, the enhanced backscatter and color ratio 
may come from the swelling of hygroscopic dust and non-dust particles in the humid air near 
clouds, or even from cloud contamination. Flowever, the depolarization ratio is expected to 
decrease and not increase near clouds, as hydrated particles tend to be more spherical than dry 
particles. Thus the apparent increase in depolarization ratio near clouds for dust volumes 
below 2 km is somewhat counter-intuitive. A possible explanation is that the hydrated and 
more spherical dust particles or those heavily mixed with marine aerosols are (mis)classified 
as non-dust aerosols due to their reduced depolarization ratio; the remaining particle 
populations will be dominated by hydrophobic dust particles that have irregular shapes and 
hence higher depolarization ratios. 

In principle, multiple scattering by undetected cloud fragments could also increase 
depolarization, but this is likely insignificant, for two reasons: (/') the increase in backscatter is 
too small to suggest strong multiple scattering near clouds, and (//') the depolarization ratio of 
all aerosols (orange curve in Figure 7d) increases only slightly near clouds, which also 
suggest that cloud contamination should be less significant if any. Another possibility could 
be that dense dust was misclassified as cloud and the observed trend could come from 
changes near thick dust, as opposed to near clouds. However, this would increase the fraction 
of dust profiles near clouds, whereas Figure 7a shows a decrease: If much of the detected 
clouds were in fact pockets of dense dust, the fraction of dust profiles should increase near 
them, as dilute dust profiles are more frequent near dense dust than far from it. 

4 Summary 

This paper uses CALIOP lidar data to examine the bulk optical properties of dust layers as 
Saharan dust moves westward over the Atlantic Ocean. It analyzes dust layers in three regions 
along the dust transport route, and examines the relationships between dust properties and the 
amount and proximity of nearby clouds. 

The study finds that the observed properties of dusty volumes are related not only to the 
meteorological conditions in the three regions, but also to the speed and duration of dry and 
wet sedimentation processes. The study examines four characteristics of dust layers: (/) the 
volume of air containing dust, (//') lidar backscatter (related to optical thickness), (//'/) color 
ratio (roughly proportional to particle size) at least for spherical particles, and (z'v) 
depolarization ratio (characterizing particle shape, with larger values for irregular particles 
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than for spherical ones). The results show that lidar backscatter and color ratio are smaller, 
while the depolarization ratio is larger in the warmer and dryer East region. 

The analysis reveals that the medians of depolarization ratio and color ratio generally increase 
with altitude in the SAL. The rate of vertical increase in depolarization ratio is significantly 
larger farther away from Africa’s west coast. 

We find the dusty volume optical properties related to cloud coverage, with backscatter and 
color ratio increase with cloudiness of surrounding areas. The effects of cloudiness are most 
prominent for dust below the SAL. The results highlight that sensitivity to cloudiness is very 
different below and within SAL. 

The results also reveal other differences between dust volume near-cloud behaviors inside and 
below the SAL. In the SAL, the fraction of aerosol samples that contain dust doesn’t depend 
on the distance to clouds, neither the median lidar backscatter, color ratio, and depolarization 
ratio. Below the SAL, the fraction of aerosol samples containing dust decreases near clouds, 
while the optical properties show noteworthy increases near clouds. The unique features of 
dust below the SAL indicate that in humid air near clouds only some large dust particles with 
much higher depolarization ratio are identified as dust by the CALIPSO detection algorithm, 
and these particles become less frequent near clouds. 

Earlier studies demonstrated that the degree of irregularity of dust affects dust optical 
properties and radiative forcing. Our observations further underline the need for assessing the 
effects of vertical separation in dust depolarization, caused by shape-dependent fall velocity, 
both in transport modeling and in estimating dust radiative forcing. In addition, our 
observations of near-cloud behaviors reveal the complexity of dust mixing with other water- 
soluble aerosols especially in the MBL, and support the hypothesis that dust, or the MBL part 
of it, becomes hygroscopic through interactions with atmospheric components in moist air, 
and this significantly affects dust optical properties. 
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4 Figure 1. Location of the 3 North- Atlantic regions examined in this study. Colors indicate the 

5 nighttime CALIOP Aerosol Optical Depth (AOD) contributed from all aerosols (Left), and 

6 contributed from dust aerosols (Right) over oceans for the June 7-July 7, 2007 period, at a 

7 pixel resolution of 2° X 2°. 
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Figure 2. Vertical profiles of (a) number of 5 km-resolution dust samples, normalized by the 
total number of dust samples at each corresponding region. The total number of dust samples 
is 45714 for E, 73141 for M and 32785 for W. In addition, the E, M, and W regions contain 
24701, 34019, and 23786 5 km-resolution vertical profiles over ocean, respectively, (b) 
median attenuated total backscatter coefficient at 532 nm, P' 532 , (c) median attenuated color 
ratio, x', (d) median volume depolarization ratio, 5'. The colors identify the profiles for each 
study region (East, Middle, and West). The error bars indicate the uncertainty of median 
values, estimated using the bootstrap algorithm. 
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3 Figure 3. Schematic illustration of Cloud Fraction (CF) definitions for 5 km resolution dust 

4 pixels. CF is the fraction of cloudy 0.3 km-resolution pixels in 5 km size areas containing 

5 dust. 
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2 Figure 4. Vertical profiles of dust properties for various cloud fractions at the 3 regions. Rows 

3 1, 2, 3 and 4 are for frequency of occurrence, attenuated total backscatter coefficient at 532 

4 nm, P'532, attenuated color ratio, and volume depolarization ratio, 5', respectively. The left, 

5 center, and right columns show the West, Middle and East regions, respectively. Results for 

6 different cloud fractions are indicated by different colors. 
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2 Figure 5. Median depolarization ratio, 5', as a function of color ratio, x' (a), and backscatter 

3 coefficient, P '532 (b). As in Figure 2, colors identify the examined regions E (black) and M 

4 (red), (c)-(d) are the same as (a) and (b), but for a different dataset from May, 25 to June, 25 

5 of 2008. 
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2 Figure 6. Rate of vertical increase in dust depolarization ratio between 3 km and 4 km 

3 altitudes, vs. distance from the African coast, represented by longitude. Particle shape- 

4 dependent differences in fall speed cause increasingly more pronounced vertical stratification 

5 as plumes move westward. 
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Figure 7. Properties of high and low dust as a function of distance to clouds, combined for the 
three regions (W+M+E): (a) fraction of detected aerosol profdes that contain dust, (b) 
attenuated backscatter coefficient at 532nm, P' 532 , (c) attenuated color ratio, (d) 
depolarization ratio, 5'. The orange curve in (a) is the number of detected “all aerosol” 
profdes as a function of distance to clouds. It is also used as denominator in calculating the 
fraction of high and low dust profiles. Orange curves in (b)-(d) show the optical properties of 
all aerosols combined. 
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Abstract. We use CALIOP nighttime measurements of lidar 
backscatter, color and depolarization ratios, as well as partic- 
ulate retrievals during the summer of 2007 to study transat- 
lantic dust properties downwind of Saharan sources, and to 
examine the influence of nearby clouds on dust. Our analy- 
sis suggests that (1) under clear skies, while lidar backscatter 
and color ratio do not change much with altitude and longi- 
tude in the Saharan Air Layer (SAL), depolarization ratio in- 
creases with altitude and decreases westward in the SAL; (2) 
the vertical lapse rate of dust depolarization ratio, introduced 
here, increases within SAL as plumes move westward; (3) 
nearby clouds barely affect the backscatter and color ratio of 
dust volumes within SAL but not so below SAL. Moreover, 
the presence of nearby clouds tends to decrease the depolar- 
ization of dust volumes within SAL. Finally, (4) the odds of 
CALIOP finding dust below SAL next to clouds are about 
2/3 of those far away from clouds. This feature, together 
with an apparent increase in depolarization ratio near clouds, 
indicates that particles in some dust volumes loose aspheric- 
ity in the humid air near clouds, and cannot be identified by 
CALIPSO as dust. 


1 Introduction 

Atmospheric mineral dust particles have significant effects 
on the climate and the environment. Despite notable ad- 
vances in modeling and satellite and ground-based measure- 
ments, dust remains the dominant factor in the uncertainty 
of aerosol radiative forcing (IPCC, 2001, 2007). Dust emit- 
ted from dry areas of Africa is transported over the North 
Atlantic Ocean to coastal areas of America with a peak of 
deposition during summer months (e.g., Prospero and Carl- 
son, 1972; Mattsson and Nihlen, 1996; Prospero and Lamb, 
2003; Torres et al., 2002; Kaufman et ah, 2005). 

The influence of dust on a radiative budget depends on 
its ability for absorbing and scattering solar and IR radiation. 
Dust optical properties are determined by the refractive index 
(i.e., chemical composition) (e.g., Sokolik and Toon, 1999; 
Wang et ah, 2002; Lafon et ah, 2006; Kahnert et ah, 2007; 
Kandler et ah, 2007, 2009; Osborne et ah, 2008; Petzold et 
ah, 2009), size and shape of dust particles (Kalashnikova and 
Sokolik, 2002; Dubovik et ah, 2006; Nousiainen, 2009). Dust 
hygroscopicity describes the particles’ ability for taking up 
water from humid air, and is rooted in the physical-chemical 
properties of dust components. Dust is mostly composed 
of water-insoluble minerals and shows nearly complete hy- 
drophobicity or at least poor hygroscopicity (e.g. Twomey, 
1977; Li- Jones et ah, 1998; Kaaden et al, 2009; Schladitz 
et ah, 2011; Ansmann et ah, 2011). However, after being 
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lifted in the air and mixed/coated with water-soluble mate- 
rials such as sea-salt, sulfate, or nitrate by atmospheric pro- 
cessing (Levin et al., 1996; Yin et al., 2002), dust hygroscop- 
icity can increase, which can cause changes in optical proper- 
ties. In addition to mineral dust particles serving as ice nuclei 
(IN) (e.g. DeMott et al., 2003; Ansmann et al. 2008), the in- 
crease in hygroscopicity may also transform dust into effec- 
tive cloud condensation nuclei (CCN) (Johnson, 1982; Wur- 
zler et al., 2000; Sassen et al., 2003; Lohmann and Feichter, 
2005; Twohy et al., 2009a) and thereby affect the formation 
and distribution of clouds and precipitation (e.g. Kelly et al., 
2007), thus also altering the radiative impact of clouds. 

The evolving shape-dependent optical and hygroscopic 
properties of dust pose a question: how does the hygroscop- 
icity of dust affect its optical properties near clouds? This 
question is pertinent as recent studies have shown that optical 
properties of clear sky aerosols are different in the vicinity 
of clouds from those far away from clouds (e.g., Clarke et al. 
2002; Twohy, et al, 2002, 2009b; Koren et al., 2007; Su et al., 
2008; Redemann et al., 2009; Tackett and Girolamo, 2009; 
Varnai and Marshak, 2011). The answer to this question is 
likely not only to improve our understanding of dust-cloud 
interactions but also yield better estimates of direct radiative 
forcing. To that end, we present analysis of dust properties 
over the North Atlantic Ocean - including near-cloud behav- 
ior - based on Cloud- Aerosol Lidar with Orthogonal Polar- 
ization (CALIOP) data (Winker, et al. 2003). 

CALIOP is a space based lidar system onboard the Cloud 
Aerosol Lidar and Infrared Pathfinder Satellite Observations 
(CALIPSO) satellite launched in 2006. CALIOP data offers 
many advantages for this study. First, since CALIOP uses 
a laser with a small footprint (~ 90 m in diameter on the 
ground), its aerosol data is not affected by the 3-D radia- 
tive enhancements of nearby clouds, which cause compli- 
cations for instruments observing reflected sunlight (Wen et 
al., 2007; Marshak et al., 2008; Varnai and Marshak, 2011). 
Second, CALIOP provides backscatter depolarization infor- 
mation at 532 nm, which allows one to distinguish (typi- 
cally) non-spherical dust from (typically) spherical droplets 
(Sassen, 2000; Murayama et al., 2001; Vaughan et al., 2004). 
Third, CALlOP’s high spatial resolution (30 m vertically and 
333 m horizontally) is well-suited for studying cloud-dust in- 
teractions that have typical scales of several kilometers (e.g., 
Koren et al., 2007; Varnai and Marshak, 2011). One concern, 
on the other hand, is that off-track clouds not detected by the 
lidar beam also affect the observed clear-sky volumes. How- 
ever, the recent study of co-located MODIS and CALIPSO 
observations demonstrates that off-track clouds do not qual- 
itatively change the aerosol property trends obtained from 
CALIOP data alone, and that the contributions from off-track 
clouds can be taken into account by using a scaling factor 
(Varnai and Marshak, 2012.) 


2 Data and methodology 

This study uses a month-long (7 June-7 July 2007) dataset 
of nighttime CALIOP Version 3 data over the North Atlantic 
Ocean (0-45° N, 0-90° W). We note that since CALIPSO or- 
bits are repeated in a 16-day cycle, our month-long dataset 
covers almost two orbital cycles, with the longitudes of clos- 
est orbits being 1.55° apart. Since in summer there are usu- 
ally multiple outbreaks per month (Huang et al., 2010), this 
month-long dataset is sufficiently long to observe the basic 
features of dust outbreaks and also sufficiently short to re- 
duce the impact of temporal changes in dust origin and mix- 
ing with other aerosols. 

CALIOP measures the total backscatter of its laser pulses 
at 532 nm and 1064 nm wavelengths, and the perpendicularly 
polarized backscatter at 532 nm. It also reports two compli- 
mentary sets of properties for each identified aerosol layer at 
5 km horizontal resolution: measured properties and derived 
particulate properties. Layer measured properties include in- 
tegrated attenuated backscatter y', integrated attenuated to- 
tal color ratio / integrated volume depolarization ratio S', 
while the corresponding layer derived particulate proper- 
ties include optical depth (r), integrated particulate color ra- 
tio (x), integrated particulate depolarization ratio (5). These 
products allow us to obtain attenuated volume backscatter 
coefficient values /i' by dividing the integrated attenuated 
backscatter values by the layer geometrical thickness. 

This paper analyzes dust observations mostly using attenu- 
ated volume properties. These properties are direct measure- 
ments without potential artifacts induced by assumptions and 
are thus suitable for describing the direct observations. On 
the other hand, since the retrieved dust particulate properties 
are free of molecular Rayleigh scattering and gas attenuation, 
they will also be used for better understanding the physics be- 
hind the observations. The results based on particulate prop- 
erties are shown in the Appendix. 

In distinguishing dust from clouds and other aerosols we 
rely on the operational CALIOP aerosol product. The op- 
erational algorithm identifies cloud and dust layers in three 
steps. 

First, it identifies particle layers based on the observed 
532 nm backscatter values (Winker et al., 2009; Vaughan et 
al., 2009). 

Second, it determines whether a detected layer is a cloud 
or aerosol layer based on its latitude, altitude, 532 nm 
backscatter, color ratio and depolarization ratio (Liu et al., 
2004, 2009). The most obvious clouds are identified at 333 m 
resolution, while the more ambiguous cases are decided at a 
coarser (1 or 5 km) resolution. Our study considers a location 
cloudy if the 1 km resolution Level 2 CALIOP cloud product 
indicates the presence of clouds. In order to reduce the im- 
pact of misclassifications between clouds and aerosols, this 
paper examines aerosol layers only if the Cloud- Aerosol Dis- 
crimination (CAD) product (Liu et al., 2004, 2009) - based 
on probability distribution functions obtained from expert 
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Fig. 1. Location of the 3 North-Atlantic regions examined in this study. Colors indicate the nighttime CALIOP Aerosol Optical Depth (AOD) 
contributed from all aerosols (left), and fraction percentage of the dust contribution to AOD (right) over oceans for the 7 June-7 July 2007 
period, at a pixel resolution of 2° x 2°. 


classifications for sample orbits - indicates that the likeli- 
hood of misclassification is less than 15%. We note, how- 
ever, that this is only the worst case, as for most layers 
used in our study the operational algorithm indicates a much 
lower probability of misclassification. Since altitude and de- 
polarization ratio are also considered in the current (Version 
3) CAD algorithm, the operational processing separates ice 
clouds from dust in this step. We note that limiting our anal- 
ysis to altitudes below 5 km and summer is also likely to re- 
duce the misclassification rate. As additional precaution, the 
statistical analysis in this paper examines median values in- 
stead of mean values, because medians are less sensitive to 
any outlying data points influenced by undetected cloud par- 
ticles. The uncertainty of median values is estimated using 
the bootstrapping algorithm (Efron and Gong, 1983). 

Third, the operational CALIOP data processing identifies 
dust layers as aerosol layers with high depolarization ratio 
values (Omar et al, 2009). Because depolarization depends 
on particle shape, it is well suited for separating typically 
non-spherical dust particles (<5 > 0.2) from usually spheri- 
cal non-dust aerosols (5 <0.075) over ocean. We note that 
this paper considers dust-containing aerosol layers identi- 
fied as either “dust” or “polluted dust” in the 5 km-resolution 
Level 2 aerosol product. (Polluted dust is dust mixed with 
biomass burning aerosols or polluted marine aerosols, with 
a depolarization ratio between those of dust and non-dust 
aerosols, 0.075 and 0.2 as discussed in Omar et al., 2009). 

To discern changes in dust properties during transatlantic 
transport, we examine dust behavior in the three regions 
shown in Fig. 1: east (E) (0-30° W), middle (M) (30-60° W) 
and west (W) (60-90° W). These three regions lie at differ- 
ent distances from the African dust sources, and cover most 
of the dust paths from Africa to America during the summer 
of 2007. 


3 Results and discussion 

The spatial and optical characteristics of African dust vary 
during the transatlantic journey (e.g., Liu et al., 2008; Huang 
et al., 2010). In this section we examine the variations in 
three steps. First, the overall statistics of dust properties in 
the three geographic regions are compared. This part focuses 
on the vertical distribution of dust samples in the CALIOP 
5 km resolution aerosol product, and on the vertical distribu- 
tions of attenuated backscatter coefficient, color ratio and de- 
polarization ratio. We then analyze the relationships between 
dust properties and cloud coverage in the three regions. Fi- 
nally, we discuss the systematic changes in dust properties 
that occur near clouds. 

In recent years there has been a large amount of research 
conducted on characterizing the physical, chemical prop- 
erties of Saharan mineral dust during its transport, such 
as Saharan Dust Experiment (SHADE) (e.g. Tanre et al., 
2003), Puerto Rico Dust Experiment (PRIDE) (e.g. Reid and 
Maring, 2003), African Monsoon Multidisciplinary Analysis 
(AMMA) (e.g. Redelsperger et al., 2006), and Saharan Min- 
eral Dust Experiment (SAMUM) (e.g. Heintzenberg et al., 
2009; Ansmann et al, 2011), etc. By using the capabilities of 
various lidars, dust profiling (e.g. Freudenthaler et al, 2009; 
Ansmann et al, 2009; Tesche et al., 2009, 2011) has provided 
insights into the vertical distributions of size, shape, extinc- 
tion and scattering of mineral dust at various sites and sea- 
sons. 

Compared to the above mentioned field experiments, this 
paper examines dust properties and the influence of clouds 
during transport based on a larger datasets that includes three 
consecutive regions along a typical transport path during 
summer. Also, the statistics is conditioned upon the presence 
of dust layers identified by CALIOP. 
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Normalized number of samples p ' 532 (km ' 1 sr' 1 ) 



Fig. 2. Vertical profiles of (a) number of 5 km-resolution dust samples, normalized by the total number of dust samples at each corresponding 
region. The total number of dust samples is 45 714 for E, 73 141 for M and 32 785 for W. In addition, the E, M, and W regions contain 24 701, 
34 019, and 23 786 5 km-resolution vertical profiles over ocean, respectively, (b) median attenuated total backscatter coefficient at 532 nm, 
^532’ me di an attenuated color ratio, x' (d) median volume depolarization ratio, S' . The colors identify the profiles for each study region 
(East, Middle, and West). The error bars indicate the uncertainty of median values, estimated using the bootstrap algorithm. 


3.1 Dust properties in the three regions 

3.1.1 Vertical distribution of dust 

Figure 2a shows the vertical distribution of dust samples in 
our three regions, normalized by the total number of dust 
samples within each region. The number of dust samples is 
defined as the number of 5 km long 270 m high volumes that, 
according to the CALIOP aerosol product, contain dust and 
have a CAD value between —70 and — 100. In the Eastern (E) 
region, more than 80 % of dust is between 1.5 km and 5.5 km 
altitude, with the peak probability around 3.5 km. This el- 
evated dust distribution is a typical result of two confining 
inversions below and above the SAL (Carlson and Prospero, 


1972). The dust remains elevated in the middle (M) region 
as well, although the mean elevation descends about 0.5 km. 
If a 3-day average transport time from Region E to M is as- 
sumed, the descending velocity from center of E to center of 
M is estimated around 1.7 mms" 1 , which is consistent with 
the typical SAL average descending velocity of 1-2 mms -1 
(Carlson and Prospero, 1972). Finally, in the West (W) region 
the chances of finding dust decrease steadily with altitude, 
and dust is rarely found above 5 km. This dramatic change in 
the vertical distribution implies that the meteorological con- 
ditions in region W are different from those sustaining the 
elevated profiles in regions E and M. Moreover, the vertical 
distribution of dust in region W indicates that (dry and wet) 
dust sedimentation has the strongest impact over region W. 
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3.1.2 Attenuated backscatter coefficient, color ratio, 
and depolarization ratio of dust 

This section examines dust optical properties in the three re- 
gions. Since dust above 5 km is rare, the analysis of dust 
optical properties will be limited to altitudes below 5 km. 
The analysis uses the Level 2 CALIOP aerosol layer product, 
which provides averaged volume properties of f$', and 5' 
as well as their particulate counterparts x and S for each dust 
layer. These layer-average values are assigned to all altitude 
bins within a dust layer when creating Fig. 2 (and 4) as well 
as Figs. Al (and A2) in the Appendix. 

The results in Figs. 2 and Al show that dust properties 
vary with altitude differently within SAL (from 1.5 to 5 km 
in altitude of regions E and M), below the SAL (below 1 .5 km 
in regions E and M), and in region W. 

In the SAL, median /J™ values are nearly constant with 
altitude, but the medians of x' and ^ as well as x and S 
increase with altitude. The increase of x' and X implies in- 
crease in size if particles are spherical. However, because of 
the complex relationship between size and color ratio for as- 
pherical particles (Bi et al., 2009), the slight increase of color 
ratio here does not imply increase of dust size. Moreover, 
further analysis in Sect. 3.2 shows the dust color ratio un- 
der clear skies to vary insignificantly with altitude. The in- 
crease of S' and S, on the other hand, suggests that aspheric- 
ity of dust increases with altitude and/or the concentration 
of mixed marine aerosols in dust volumes decreases with al- 
titude, since large asphericity of dust particles implies large 
S. 

Below the SAL in the regions E and M, the median S' (as 
well as S) increases with altitude. This is the result of dust 
mixing with non-dust marine aerosols in the moist air con- 
fined between the marine surface and the inversion created by 
the dry and warm SAL aloft. Since the concentration of wet 
marine aerosols is much higher below than inside the SAL, 
backscatter from these aerosols contributes significantly to 
the lidar signals and reduce the depolarization ratios of dust 
volumes below the SAL. 

In region W, the median of /L 32 decreases with altitude 
and the medians of x' and ( as well as x and <5) increase 
with altitude. Similar to the case in the E and M, the increase 
in color ratio does not necessarily imply the increase of dust 
size, while the increase of depolarization ratio indicates the 
increase of dust asphericity and/or decrease of concentration 
of marine aerosols mixed in dust volumes at higher altitude. 
The median of /3/-p decreasing with altitude is the direct re- 
sult of reduced dust and/or marine aerosol concentration at 
higher altitudes. 

A comparison of dust properties in regions E, M, and W 
also reveals several features at different transport stages. For 
example, at most altitudes the medians of and y' in- 
crease westward, whereas the median of S' (and of x and S ) 
tends to decrease westward. We note that the further study 
of fully cloud-free regions in Sect. 3.2 indicates westward 


/l™ increases only below ~ 2 km, with no significant west- 
ward changes above 2 km. The westward increase of /3L 2 
below 2 km is likely the result of mixing with higher concen- 
tration of marine aerosols in the W than in E. Monthly data 
from the Goddard Earth Observing System Model Version 5 
(GEOS-5) Modern Era Retrospective-Analysis for Research 
and Applications (MERRA) for June 2007 shows increasing 
surface wind speed westward along the dust path in Figure 
1, with values around ~13ms _1 , ~9ms -1 and ~2ms _1 
at 850.0 hPa at 15° N latitude for the W, M, and E regions. 
A westward increase in wind speed is also observed by Da- 
vidi et al. (2012). The stronger surface wind speed generates 
higher concentration of marine aerosols and better mixing 
with dust towards the W. This is consistent with the observed 
westward trend of decreasing particulate 8. 

The westward decreasing trend of depolarization ratio be- 
low the SAL also suggests that the amount of marine aerosols 
mixed with dust in the Marine Boundary Layer (MBL) in- 
creases westward. This issue will be further explored in 
Sect. 3.2, which also considers the influence of cloudiness. 

However, some features in Fig. 2 cannot be explained by 
contributions from non-dust marine aerosols in dust sample 
volumes, and are likely caused by changes in the properties 
of dust particles instead. For example, Fig. 2d shows that 
above 3.5 km, the median S' is larger in region M than in 
region E. The features remain similar for S in as shown in 
Fig. Al. This cannot be explained by mixing from below or 
more well-mixing westward, because the mixing at high al- 
titude of SAL (above 3.5 km) is rare. Instead, the observed 
tendency is likely related to lower fall speed for aspherical 
dust particles: as the more spherical particles fall faster, this 
leaves an increasingly non-spherical dust population at high 
altitudes as the air moves to region M. The plausibility of this 
scenario is also supported by simulations for highly irregu- 
lar particles falling slower than more spherical ones, because 
of greater air resistance (Ginoux, 2003). This issue will be 
further discussed in Sect. 3.3.2. 

We note that the slight increase of x' or decrease of x 
of dust volumes from E to M (in Figs. 2c and Al) appears 
different from the behavior of dust Angstrom exponents re- 
trieved by MISR and MODIS, the latter displaying no sig- 
nificant changes during transatlantic transport (Kalashnikova 
and Kahn, 2008). This apparent difference can be removed 
when we examine the color ratios under clear skies as dis- 
cussed in Sect. 3.2, in which no significant difference in color 
ratio is found in all 3 regions. 

3.2 Correlation of dust properties with cloud fraction 

This section examines the relationships between dust prop- 
erties and cloudiness in the three study regions. We charac- 
terize cloudiness through the cloud fraction (CF), defined for 
each dust-containing 5 km-size column as the ratio of number 
of cloudy 0.333 km profiles to the total number of 0.333 km 
profiles in the column. Simply put, if the number of cloudy 
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0. 333 km profiles is m, the cloud fraction is m 15 _1 . We note 
that in addition to the cloud fraction varying between 0 and 

1, the relative location of dust and clouds within 5 km wide 
columns can also vary (Fig. 3). We also note that unlike the 
conventional cloud fraction that is based on 2-dimensional 
(2-D) images, our definition here is based on 1-dimensional 
(1-D) measurements along the CALIPSO track. Although 
off-track clouds may influence dust properties along the track 
even for CFi_d = 0, there is direct statistical relationship be- 
tween the defined CFi_d here and the conventional CF 2 -D 
(Varnai and Marshak, 2012), and CFi_d is still a generally 
useful indicator of cloud coverage. 

The results in Fig. 4 show that dust properties are closely 
related to CF in all three regions. The main features of the 
relationship are as follows. 

First, the top row of Fig. 4 reveals that a smaller fraction 
of dust samples occurs under clear skies in region M than 
in region E. This is because the SAL is warmer and drier in 
the East, and so the conditions are less favorable for cloud 
formation in region E than region M. 

Second, rows 2 and 3 in Fig. 4 reveal that within each re- 
gion, the median values of f }' 5 3? and x' are larger for higher 
CF’s. This feature is likely caused by aerosols getting hy- 
drated and swelling in humid regions containing clouds, al- 
though undetected cloud particles may also contribute. The 
figure also shows that in regions E and M, the increase in 
backscatter and color ratio is more pronounced below the 
SAL than inside it. The swelling is greater below the SAL 
than inside it both because clouds and high humidity are 
more common below the SAL, and because hygroscopic ma- 
rine aerosols are fairly abundant at low altitudes even in dust 
layers, whereas the SAL is dominated by less hygroscopic 
dust particles. 

Third, within each region, two opposite trends of corre- 
lations appear between &' and CF : inside the SAL, the me- 
dian <5' of dust is always larger in clear sky than in cloudy 
skies; whereas below the SAL, the median S ' of dust is al- 
ways smaller in clear sky than in cloudy skies. The domains 
of these opposite behaviors can be separated in the fourth 
row of Fig. 4 roughly at the crossing point of the red curve 
(CF = 0) and the blue curve (0 < CF < 0.6). We note that 
these crossing points are approximately at the altitude of the 
bottom of the SAL. The opposite trends inside and below the 
SAL clearly indicate a different dust volume depolarization 
response to increased humidity. The possible mechanisms af- 
fecting the apparent depolarization ratio of dust volumes be- 
low the SAL will be discussed in Sect. 3.4. 

It should be also noted that since the dust layers discussed 
in this subsection include those below clouds, the quality 
of lidar signals returned from such dust may have been de- 
graded by the clouds above. This can further affect the re- 
trieved dust particulate properties seen in Fig. A2. Therefore 
the discussion of CF influence is based on using the volume 
property here instead of using retrieved particulate property. 
However, since signals from dust layers under clear skies are 
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Fig. 3. Schematic illustration of Cloud Fraction (CF) definitions for 
5 km resolution dust pixels. CF is the fraction of cloudy 0.3 km- 
resolution pixels in 5 km size areas containing dust. 


not affected by clouds above, particulate properties (Fig. A2) 
for CF = 0 reflect dust properties without quality degrading 
effects from clouds. 

3.3 Features of dust volumes in the SAL under 
clear skies 

Unlike the dust volumes below the SAL, where the dust is 
mixed with humidified non-dust aerosols, the SAL is domi- 
nated by dust particles. This subsection examines several fea- 
tures of dust volumes inside the SAL in regions E and M. To 
reduce the effects of clouds, dust volumes are limited to only 
those under clear skies. In addition, region W is excluded 
because of its low number of dust samples inside the SAL 
(Fig. 2a). 

3.3.1 Relationship of depolarization ratio with 
color ratio 

Figures 2 and 4 show that while backscatter is fairly uniform 
vertically, both color ratio and depolarization ratio increase 
markedly with altitude inside the SAL of regions E and M. 
These coinciding increases suggest systematic relationship 
between the depolarization and color ratio (Fig. 5). Figure 5a 
reveals a positive relationship between the depolarization ra- 
tio and color ratio of dust volumes in the SAL of regions E 
and M. In addition, the dynamic ranges of depolarization ra- 
tio and color ratio are much wider in region M than in region 
E. Figure 5b confirms that similar relationships are valid for 
a different dataset (covering 25 May-25 June 2008) as well. 

The relationships shown in Fig. 5 and the similarity of 
results from the two independent datasets can be attributed 
to the steady altitude-dependence of color ratio and depo- 
larization ratio inside the SAL. As discussed in Sect. 3.1.2, 
these altitude dependences are likely caused by two mech- 
anisms: (i) a decrease with altitude in the concentration of 
non-dust particles mixed in from below, and (ii) different fall 
speeds vertically separating the relatively more spherical dust 
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Fig. 4. Vertical profiles of dust properties for various cloud fractions at the 3 regions. Rows 1, 2, 3 and 4 are for frequency of occurrence, 
attenuated total backscatter coefficient at 532 nm, /Sj 3? , attenuated color ratio, and volume depolarization ratio, S', respectively. The left, 
center, and right columns show the West, Middle and East regions, respectively. Results for different cloud fractions are indicated by different 
colors. 
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Fig. 5. Median depolarization ratio, S', as a function of color ratio, x' ( a )- As in Fig. 2, colors identify the examined regions E (black) and M 
(red), (b) is the same as (a), but for a different dataset from 25 May to 25 June 2008. 


particles from the least spherical ones. This latter mechanism 
dominates and is further explored in Sect. 3.3.2. In addition, 
we note that the depolarization ratio is predominantly influ- 
enced by asphericity of the dust particles (e.g. Ansmann et 
al., 2003). 

3.3.2 Relationship between the vertical increase in 
depolarization ratio and longitude 

As indicated in Figs. 2 and Al, the depolarization ratio not 
only increases with altitude in the SAL, but also has a larger 
increase rate in region M than E. As mentioned above, the in- 
crease may come from more spherical and less spherical dust 
particles getting vertically separated because of their differ- 
ent sedimentation speeds. The upward increase in S' and 8 
could then be stronger in region M simply because the sedi- 
mentation process has more time to work by the time the dust 
reaches region M. 

Numerous studies have demonstrated that the fall speed of 
atmospheric particles depends on their shape (e.g., Cheng et 
al., 1988; Ginoux, 2003) as well as size. Since particles with 
irregular shapes have greater cross-sectional areas and drag- 
coefficients, they experience stronger drag force in the air - 
which implies that more irregular particles of the same mass 
fall slower than more spherical ones. Note that a sphere is 
the most compact object (least surface area for a given vol- 
ume) and it experiences least drag for a given mass. (Here 
we assume that dust particle shape does not change sys- 
tematically with particle size, which also greatly impacts 
fall speed. This assumption is plausible because field exper- 
iments show that distribution of aspect ratio is weakly re- 
lated to particle size; e.g. Chou et al., 2008; Handler et al., 
2009.) As a result, shape-induced vertical separation will en- 
sue as dust is advected westward, with irregular particles in- 
creasingly predominant in the upper portions of SAL. At a 


constant altitude, this stratification is expected to widen the 
dynamic range of depolarization ratios with downstream dis- 
tance from the dust source. This is indeed the case. 

To that end, we divide regions E and M into sub-regions 
covering 10° wide longitude bands, and examine the average 
difference between the depolarization ratios at 3 and 4 km al- 
titudes for each region. As shown in Fig. 6 (and Fig. A4), 
the average difference S' Akm - S ' 3 ktn (or S 4 km - <$ 3 km) keeps 
increasing with the distance from the west coast of Africa. 
This result implies that the observed change in volume de- 
polarization ratio within SAL is most likely caused by the 
greater drag of aspherical dust particles. In addition, this re- 
sult is consistent with the observation of higher aspect ratio 
of long rang transported Saharan dust Reid et al. (2003). 

3.4 Dust volume properties near clouds 

Relative humidity usually increases as clouds are approached 
and this causes nearby aerosols to swell and get hydrated 
(acquire thin film of water) or even activated as haze (e.g., 
Twohy et al., 2009b). Observing changes of dust character- 
istics near clouds can help improve our understanding of the 
effect of high relative humidity and clouds on dust particles. 
Figure 4 has shown that the backscatter /T and color ratio 
x' of dust volumes increase with cloud fraction both in and 
below the SAL, whereas the depolarization ratio S' changes 
with cloud fraction differently in and below the SAL. This 
finding indicates that dust properties in and below the SAL 
are different. This section further examines the near cloud 
behaviors of dust volumes in and below the SAL. We note 
that although the base altitude of the SAL may vary during 
westward transport, as shown in Figs. 2 and 4, this analy- 
sis uses constant separation altitude of 2 km for convenience. 
Since there is no aerosol particulate retrieval available at res- 
olutions higher than 5 km, our analysis uses CALIOP Level 1 
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Fig. 6. Rate of vertical increase in dust depolarization ratio between 
3 and 4 km altitudes, vs. distance from the African coast, repre- 
sented by longitude. Particle shape-dependent differences in fall 
speed cause increasingly more pronounced vertical stratification as 
plumes move westward. 


data to examine changes in backscatter, color ratio, and de- 
polarization as a function of distance to clouds at a resolution 
of 0.333 km. In this analysis a 0.333 km resolution clear sky 
profile is considered a dust profile if it is included in one or 
multiple 5 km resolution dust layer(s). Moreover, a clear sky 
requirement of no 0.333 km, 1.0 and 5 km clouds above or 
below used in selecting dust profiles near clouds avoids the 
signal degradation mentioned in the previous subsection ear- 
lier. 

Figure 7 illustrates the behavior of dust as a function 
of distance to clouds. The orange curve corresponds to all 
aerosol samples while the black and green ones to high (in 
the SAL) and low (below the SAL) dust, respectively. Fig- 
ure 7a shows that the fraction of dust profiles over all de- 
tected aerosol profiles decreases dramatically near clouds for 
dust at altitudes below the SAL, but remains relatively stable 
for those in the SAL. The stable behavior in the SAL can be 
explained by the fact that low clouds confined to the bound- 
ary layer by the inversion at the base of the SAL have little 
impact on humidity inside the SAL. The near-cloud drop in 
the fraction of dust profiles below the SAL may result from 
several factors. First, the chances of wet removal are higher 
near clouds, and this can lower the fraction of dust profiles. 
Second, swelling in the humid air near clouds makes parti- 
cles more spherical (especially if water-soluble particles pol- 
lute dust crystals), resulting in the reduction of depolariza- 
tion ratio thus the dust signature (Omar et al., 2009); con- 
sequently, some dust-containing profiles are (mis)classified 
as non-dust aerosol. Assuming that below 2 km altitude, the 


fraction of dust profiles is constant beyond 5 km from clouds 
(as most humidity changes occur within 5 km from clouds), 
the roughly 2/3 drop in the fraction of dust profiles near 
clouds implies that at least 2/3 of dust profiles in the MBL 
are polluted and hygroscopic. This results in CALIOP miss- 
ing dust in about 1 /3 of dust profiles that occur within 5 km 
from clouds. 

Figure 7b an d also show that backscatter, color ratio, and 
depolarization ratio all increase near clouds for dust layers 
below 2 km, but they remain fairly constant for dust lay- 
ers above 2 km. The stable behavior in the SAL occurs be- 
cause most clouds are below the SAL and have little impact 
on dust in the SAL. In addition, the dust population in the 
SAL is dominated by hydrophobic particles. For dust vol- 
umes below 2 km, the enhanced backscatter and color ratio 
may come from the swelling of hygroscopic dust and non- 
dust particles in the humid air near clouds, or even from 
cloud contamination. However, the depolarization ratio is ex- 
pected to decrease and not increase near clouds, as hydrated 
particles tend to be more spherical than dry particles. Thus 
the apparent increase in depolarization ratio near clouds for 
dust volumes below 2 km is somewhat counter-intuitive. This 
counter-intuitive feature of 8 ' near clouds remains the same 
for the estimated particulate depolarization ratio 8 in Fig. A5. 
A possible explanation is that the hydrated and more spheri- 
cal dust particles or those heavily mixed with marine aerosols 
are (mis)classified as non-dust aerosols due to their reduced 
depolarization ratio; the remaining particle populations will 
be dominated by hydrophobic dust particles that have irregu- 
lar shapes and hence higher depolarization ratios. 

In principle, multiple scattering by undetected cloud frag- 
ments could also increase depolarization, but this is likely 
insignificant, for two reasons: (i) the increase in backscatter 
is too small to suggest strong multiple scattering near clouds, 
and (ii) the depolarization ratio of all aerosols (orange curve 
in Fig. 7d) increases only slightly near clouds, which also 
suggest that cloud contamination should be less significant 
if any. Another possibility could be that dense dust was mis- 
classified as cloud and the observed trend could come from 
changes near thick dust, as opposed to near clouds. However, 
this would increase the fraction of dust profiles near clouds, 
whereas Fig. 7a shows a decrease: if much of the detected 
clouds were in fact pockets of dense dust, the fraction of dust 
profiles should increase near them, as dilute dust profiles are 
more frequent near dense dust than far from it. 


4 Summary 

This paper uses CALIOP lidar data to examine the bulk opti- 
cal properties of dust layers as Saharan dust moves westward 
over the Atlantic Ocean. It analyzes dust layers in three re- 
gions along the dust transport route, and examines the rela- 
tionships between dust properties and the amount and prox- 
imity of nearby clouds. 
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Fig. 7. Properties of high and low dust as a function of distance to clouds, combined for the three regions (W + M + E): (a) fraction of detected 
aerosol profiles that contain dust, (b) attenuated backscatter coefficient at 532 nm, /I 539 , (c) attenuated color ratio, x' , (d) depolarization ratio, 
S'. The orange curve in (a) is the number of detected “all aerosol” profiles as a function of distance to clouds. It is also used as denominator 
in calculating the fraction of high and low dust profiles. Orange curves in (b-d) show the optical properties of all aerosols combined. 


The study finds that the observed properties of dust vol- 
umes are related not only to the meteorological conditions in 
the three regions, but also to the speed and duration of dry 
and wet sedimentation processes. The study examines four 
characteristics of dust layers: (i) the volume of air contain- 
ing dust, (ii) lidar backscatter (related to optical thickness), 
(iii) color ratio (related to particle size at least for spher- 
ical particles), and (iv) depolarization ratio (characterizing 
particle shape, with larger values for irregular particles than 
for spherical ones). The results show that lidar backscatter 
and color ratio under clear skies do not change much dur- 
ing transport, while the depolarization ratio is larger in the 
warmer and dryer East region. 

The analysis reveals that the medians of depolarization ra- 
tio generally increase with altitude in the SAL. The rate of 
vertical increase in depolarization ratio is significantly larger 
farther away from Africa’s west coast. 

We find that the optical properties of dust volumes are re- 
lated to cloud coverage, with backscatter and color ratio in- 


creasing with the cloudiness of surrounding areas. The ef- 
fects of cloudiness are most prominent for dust below the 
SAL. The results highlight that sensitivity to cloudiness is 
very different below and within SAL. 

The results also reveal other differences between dust vol- 
ume near-cloud behaviors inside and below the SAL. In the 
SAL, the fraction of aerosol samples that contain dust doesn’t 
depend on the distance to clouds, the median lidar backscat- 
ter, color ratio, and depolarization ratio. Below the SAL, 
the fraction of aerosol samples containing dust decreases 
near clouds, while the optical properties show noteworthy 
increases near clouds. The unique features of dust below the 
SAL indicate that in humid air near clouds only some large 
dust particles with much higher depolarization ratio are iden- 
tified as dust by the CALIPSO detection algorithm, and these 
particles become less frequent near clouds. 

Earlier studies demonstrated that the degree of irregularity 
of dust affects dust optical properties and radiative forcing. 
Our observations further underline the need for assessing the 
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effects of vertical separation in dust depolarization, caused 
by shape-dependent fall velocity, both in transport modeling 
and in estimating dust radiative forcing. In addition, our ob- 
servations of near-cloud behaviors reveal the complexity of 
dust mixing with other water-soluble aerosols especially in 
the MBL, and support the hypothesis that dust, or the MBL 
part of it, becomes hygroscopic through interactions with at- 
mospheric components in moist air, and this significantly af- 
fects dust optical properties. 


Appendix A 

CALIOP measures backscattered signals at two wavelengths, 
532 nm and 1064 nm, and also measures the perpendicularly 
polarized return signal at 532 nm. For each detected aerosol 
layer, the CALIOP Level 2 aerosol product reports these 
measurements as the “volume properties” examined in the 
main body of this paper. Since, however, the directly mea- 
sured signals are affected by molecular Rayleigh scattering 
and by two-way attenuation along the path to the observed 
aerosol layers, the CALIOP Level 2 aerosol product also in- 
cludes “particulate properties”. These particulate properties 
are obtained by estimating the contributions of Rayleigh scat- 
tering and ozone absorption using GEOS-5 molecular den- 
sity profiles, and considering aerosol and cloud attenuation 
above using appropriate typical lidar ratio values. 

Since aerosol particulate properties are free of molecular 
Rayleigh scattering and gas or cloud attenuation, observing 
the variations in dust particulate properties can help under- 
stand changes in dust properties during transatlantic trans- 
port and near clouds. This appendix presents dust particulate 
color ratio and depolarization ratio results that correspond to 
the volume property results discussed in Sect. 3. For conve- 
nience, the figure numbering is parallel to those in main text. 

As expected from the volume products, particulate color 
ratios in the three regions (Fig. Al left) show increasing trend 
with altitude. However, they also show a decreasing trend 
westward that is opposite to westward increasing trend of 
volume color ratio (Fig. 2a). Particulate depolarization ratio 
values (Fig. Al right) show similar trends with altitude and 
longitude as volume depolarization ratios do (Fig. 2b). 

Figure A2 reveals that the dependence on cloud cover - 
that is, the difference between the red curves representing 
cloud free areas and the other curves - is markedly differ- 
ent for particulate properties than for the volume properties 
shown in Fig. 4. Dust particulate color ratio is larger for clear 
skies (CF = 0) than for cloudy skies (CF ^0), while volume 
color ratio is larger for cloudy skies than for clear skies. Dust 
particulate depolarization ratio is always larger for clear skies 
than for cloudy skies, while volume depolarization ratio is 
larger in clear skies than cloudy skies only above ~ 2 km. 

The particulate properties’ dependence on cloudiness may 
come from swollen non-dust aerosols or undetected cloud 


particles occurring in dust volumes in cloudy regions: these 
large spherical particles can increase color ratio and decrease 
depolarization ratio. We note, however, that if a dust layer 
is below a cloud layer, attenuation in the cloud above can 
reduce quality of particulate products. The uncertainties are 
likely larger for dust particulate properties than for volume 
properties because of uncertainties in the lidar ratios used for 
taking cloud attenuation into account when obtaining partic- 
ulate properties. Under clear skies, however, dust measure- 
ments are not affected by clouds above, and particulate prop- 
erties are likely as accurate as volume properties; their dif- 
ference comes mainly from well known Rayleigh scattering 
and ozone absorption. 

Figure A3 shows similar positive relationships between 
the particulate color ratio and depolarization ratio as in Fig. 5. 
These particulate property relationships under clear skies 
confirm the relationships expressed in the manuscript. 

Figure A4 shows a similar increasing depolarization in- 
crease rate to that in Fig. 6. This indicates that the increase 
depolarization ratio lapse rate does not come from clouds, 
Rayleigh scattering, or ozone absorption. 

Since the CALIOP Level 1 product does not provide re- 
trieved dust particulate depolarization ratio values, we exam- 
ine near-cloud changes in dust particulate depolarization ra- 
tio values calculated using the following formula (e.g., Cairo 
et al, 1999; Liu et al, 2008; Omar et al, 2009) 

s 8 Y [R + R-S m ~S m ]-S m 

Sp= s _i + s . Jm _ s , (A1) 

where <5 p , <S V and <5 m are the particulate depolarization ratio, 
volume depolarization ratio, and molecular depolarization 
ratio, while R is the backscatter ratio determined by the par- 
ticulate backscatter coefficient /3 p and molecular backscatter 
coefficient using R = (fi p + 

The estimated dust particulate depolarization ratio in 
Fig. A5 shows qualitatively the same near cloud behavior as 
Fig. 7d. This confirms that the behavior in Fig. 7d is not the 
result of Rayleigh scattering or ozone attenuation. 
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Fig. Al. Vertical distribution of retrieved dust particulate color ratio (left) and depolarization ratio (right). These plots correspond to Fig. 2b 
and d. 


w 


M 




Fig. A2. Effect of cloud fraction on the vertical distribution of dust particulate color ratio (first row) and depolarization ratio (second row). 
These two rows correspond to the 3rd and 4th rows of Fig. 4. 
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Fig. A3. Depolarization ratio as a function of color ratio, based on dust particulate properties. These panels correspond to Fig. 5a and b 
according from Left to Right. 



Longitude, deg 


Fig. A4. Rate of vertical increase in dust particulate depolarization 
ratio between 3 and 4 km altitudes, vs. distance from the African 
coast. This corresponds to the Fig. 6. 
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Fig. A5. Estimated particulate depolarization ratio as a function of 
distance to clouds for low dust. This panel corresponds to Fig. 7d. 
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